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Abstract. We observed the middle-aged radio pulsar 
B0656-I-14 with a CCD detector at the 6-m telescope. 
Broadband BVRI images show the following magnitudes 
of the pulsar counterpart: B = 24.85(4-0.19, —0.16), V = 
24.90(-b0.16,-0.14), R = 24.52(-f0.12,-0.11), I = 
23.81(-l-0.27, —0.21). We fitted the UV-optical (space -I- 
ground-based) data with a two-component model which 
combines a power law (non-thermal component) with a 
thermal spectrum emitted by the neutron star surface. 
The power law component, with the energy power-law 
index a = 1.5(-|-1.1,—1.2), dominates in the observed 
range. Constraints on the thermal component correspond 
to the Rayleigh-Jeans parameter G = T^{Rw/d^m)^ = 
4.1(4-2.1,—4.1), where T = lO^Te K is the brightness 
temperature, Rao = 10i?io km is the neutron star ra¬ 
dius as seen by a distant observer, and d = SOOdsoo pc 
is the distance. The shape of the optical-UV spectrum of 
PSR B0656-I-14 differs considerably from those observed 
from other pulsars: the middle-aged Geminga and young 
Crab, Vela. 

Key words: stars: neutron - pulsars: individual: 
PSR B0656-I-14 - techniques: photometric - techniques: 
image processing 


1. Introduction 


Optical radiation has been detected from only ~ 1% 
of ~ 750 radio pulsars (see Bignami & Caraveo 1996 
for a recent review). Optical counterparts were reported 
for pulsars of quite different ages, f rom the very young 
PSR B0531-I-21 (Crab pulsar — see Percival et al. 199'^ 
to old PSR B1929-I-10 and B0950-I-08 (Pavlov, Stringfel- 
low & Cordova 1996a). Of special interest among these 
objects are the middle-aged (r ^ 10® — 10® yr), soft X- 
ray pulsars PSR B0656-I-14, J0633-I-1746 (Geminga) and 


B1055-52 studied extensively with the X-ray observatories 
ROSAT and ASCA. Their thermal-like X-ray radiation 
is believed to be emitted by surface layers of the neutron 
stars (e. g., Ogelman 1995), so that studying this radiation 
enables one to estimate surface temperatures of neutron 
stars (NSs) of different ages, to compare them with the 
current models for the NS cooling, and to constrain prop¬ 
erties of the superdense matter in the NS interiors (e. g., 
Pavlov et al. 1995). Since the surface temperature, as well 
as other parameters of radiating layers, can be measured 
much more accurately if the surface radiation is also ob¬ 
served in the optical-UV range (Pavlov et al. 1996b), in¬ 
vestigation of the optical-UV radiation of these objects is 
potentially very important. 

In this paper we report multi-color CCD observa¬ 
tions of the middle-aged PSR B0656-I-14 (r = P/2p = 
1.1 X 10® yr) with the 6-meter telescope (BTA) at the Spe¬ 
cial Astrophysical Observatory of the Russian Academy 
of Sciences (SAO RAS). The optical counterpart of this 
pulsar was first detected with the ESO NTT and 3.6- 
meter telescopes in the V band {V ~ 25) by Caraveo et 
al. (1994a) at a 3a level. Pavlov et al. (1996a) observed 
the pulsar with the Faint Object Camera (FOC) onboard 
the Hubble Space Telescope (HST) in a longpass filter 
F130LP (AA = 2310 — 4530 A). They confirmed the opti¬ 
cal detection at a very high confidence level {S/N = 52) 
and concluded that the observed radiation is mainly of a 
nonthermal (most likely, magnetospheric) origin, although 
a contribution of a thermal (Rayleigh-Jeans) component 
emitted by the NS surface cannot be excluded. In order 
to separate the thermal component, and to understand 
the nature of the nonthermal radiation, the flux of the ob¬ 
ject should be measured in several spectral bands covering 
as broad spectral range as possible. Mignani, Caraveo & 
Bignami (1997) observed PSR B0656-I-14 with the HST 
Planetary Camera (PC) in the F555W filter and deter¬ 
mined its V magnitude more accurately: V = 25.1 4i 0.1. 
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PSR 0656+14 


V 2400s 


Fig. 1. Image of the PSR B0656+14 field in the V filter with accumulation time of 2400 sec. The center of the 30” x 30" 
fragment coincides with the pulsar position. The numbers denote the stars whose magnitudes are given in Table 3. 


Obviously, observations with other filters are needed to 
find the shape of the spectral flux of the optical counter¬ 
part. Thus, we included such observations in our general 
program of deep broad-band CCD photometry of nearby 
pulsars of the northern sky with the 6-meter telescope. 

In our first BTA observations of PSR B0656-I-14 (Jan¬ 
uary 1996) in the B, V and i? filters (Kurt et al. 1997a) we 


obtained estimates of its brightness in the V and, for the 
first time, B and i? bands and confirmed the nonthermal 
nature of its optical spectrum. However, because of poor 
weather conditions during those observations, the fluxes 
were determined with substantial uncertainties. Subse¬ 
quent observations of the same object (November 1996) 
were more successful (Kurt et al. 1997b). We describe 
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Table 1. Observations of PSR B0656+14 on November 11/12, 
1996. 


time 

exp. 

exposure 

seeing 

Z 

backgr. 

UT 

ID 

(sec) 

(arcsec) 

(deg) 

(mag/ 

arcsec^ 

22:37 

B 1 

600 

1.64 

40.63 

22.41 

22:49 

V 1 

600 

1.56 

39.00 

21.60 

23:01 

R 1 

600 

1.42 

37.43 

20.70 

23:11 

I 1 

300 

1.31 

35.95 

19.39 

23:18 

I 2 

300 

1.36 

35.17 

19.34 

23:25 

B 2 

600 

1.69 

34.42 

22.39 

23:36 

V 2 

600 

1.42 

33.21 

21.61 

23:47 

R 2 

600 

1.31 

32.13 

20.80 

23:59 

I 3 

300 

1.42 

31.22 

19.52 

00:06 

I 4 

300 

1.33 

30.76 

19.47 

00:12 

I 5 

300 

1.26 

30.39 

19.45 

00:19 

B 3 

600 

1.67 

30.06 

22.34 

00:30 

V 3 

600 

1.31 

29.65 

21.52 

00:42 

R 3 

600 

1.45 

29.43 

20.64 

01:13 

R*4 

600 

1.45 

29.88 

20.64 

01:25 

I 6 

300 

1.37 

30.42 

19.29 

01:31 

I 7 

300 

1.40 

30.50 

19.32 

01:41 

I 8 

300 

1.30 

31.57 

19.41 

01:48 

B 4 

600 

1.75 

32.10 

22.27 

01:59 

V 4 

600 

1.61 

33.18 

21.41 

02:11 

R 5 

600 

1.51 

34.39 

20.57 

02:22 

R 6 

600 

1.53 

35.74 

20.53 

02:34 

I 9 

300 

1.51 

37.20 

19.22 

02:41 

I 10 

300 

1.53 

38.05 

19.05 

02:47 

I 11 

300 

1.51 

38.95 

18.75 


these observations and their reduction in Section 2 and 
discuss the properties and origin of the optical spectrum 
of PSR B0656+14 in Section 3. 

2. Observations and data analysis 

2.1. Observations 

Photometric observations of PSR B0656+14 were carried 
out with the 6 -meter telescope at SAO RAS on 11/12 
November, 1996 with a CCD photometer installed at the 
Primary Focus. We used the “Electron ISD017A” CCD; 
its format of 1040 x 1160 pixels corresponds to the field 
of view 2'.38 x 2'.66. The CCD photometer was employed 
in the 2x2 binning mode, so that each of the 520 x 580 
zoomed pixels (referred to as ‘pixels’ hereafter) has the 
angular size of 0/274 x 0/274. We used the gain of 2.3e“ 
per DN (Data Number) pixel. The shutter timing accuracy 
is better than 0.1 sec. 

The observations were carried out with filters close to 
the B, V, R and I filters of the Cousins system. Table 1 
gives basic parameters of observational conditions: start¬ 
ing time of each exposure, filter and exposure number, 
duration of exposure, seeing, zenith distance Z, and sky 


background in stellar magnitudes per arcsec^. 


Table 2. Coordinates of objects in the PSR B0656-I-14 field. 


No. 

01950 

^1950 

relative 

coordinates 

(arcsec) 

distance 
from PSR 
(arcsec) 

1 

06:57:00.97 

14:18:51.5 



2 

06:56:57.62 

14:19:05.1 



3 

06:56:55.63 

14:17:43.6 



4 

06:57:00.02 

14:17:36.6 



5 

06:57:01.51 

14:18:24.1 



6 

06:56:56.63 

14:18:31.3 



ol,PSR 

06:56:57.84 

14:18:34.2 

206.5; 203.0 

0.0 

o2 

06:56:57.75 

14:18:36.2 

207.5; 204.3 

1.7 

o3 

06:56:57.87 

14:18:37.4 

206.1; 207.1 

4.1 

o4 

06:56:57.47 

14:18:34.7 

212.0; 204.0 

5.1 

o5 

06:56:57.58 

14:18:34.7 

209.5; 204.0 

3.1 


2.2. Data reduction 

Standard data reduction includes subtraction of the 
“bias”, — an additional component of the CCD signal, 
flat-fielding (correction for non-uniform sensibility of the 
detector elements), and removing of space particle traces. 
The data were processed with the use of the MIDAS soft¬ 
ware. Processing involved all the obtained images except 
for R*4 where the space particle trace was in the immedi¬ 
ate vicinity of the pulsar position. Figure 1 shows the field 
around the pulsar (sum of 4 exposures in the V filter). 

2.3. Astrometrical referencing 

According to Thompson and Cordova (1994), the VLA 
position of the pulsar is oigso = 06^ 56’” 57/942; ^1950 = 
14‘^18'33/80 (epoch 1992.98); the formal uncertainties of 
the radio position ((y.'09 and 0/03 for a and 6, respectively) 
do not include possible systematic errors which may be as 
large as ~ 0/2 (cf. Pavlov et al. 1996a). With allowance 
for the proper motion of the pulsar, g,a = -1-0.73 ± 20 
and fis = —26± 13 milliarcsec yr“^ (Pavlov et al. 1996a), 
the expected position of the pulsar at the epoch of our 
observations (1996.86) is 01959 = 06^ 56’” 57f96(±0/02), 
J 1950 = 14M8'33'/7(±0'/2). 

The position of the candidate optical counterpart rela¬ 
tive to nearest bright stars found with the use of the HST 
Guige Star Catalog (CSC) has been presented by Caraveo 
et. al. (1994b). The astrometric referencing of our data is 
based on coordinates taken for the field stars around the 
pulsar from the Digitized Sky Servey (nominal accuracy 
^ 1 "). 
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Table 2 gives the coordinates of the pulsar optical can¬ 
didate, 6 bright stars (denoted by the numbers 1 through 
6 in Fig. 1), and 4 objects within ~ 5" around the pulsar. 
For these objects, the coordinates corresponding to the 
frames of Figs. 3, 4 and 6 are also presented. 

2.^. Photometry 

Photometric referencing was carried out by observing four 
standard stars from the Landolt catalog (1992): 112822, 
GD71 and RU 149 (8 stars) and star 2015 from Weis 
(1991). Photometric conditions remained stable during the 
whole night, with an insignificant increase of absorption 
at the blue part of the spectrum at the end of the observa¬ 
tions. For the extinction factors, we used average values of 
the observatory: kb = 0.34, ky = 0.21, kr = 0.15, ki = 0.10 
(in stellar magnitudes). From observations of the stan¬ 
dard stars we derived the following transition equations 
from the instrument system (magnitudes b, v, r, i) to the 
Cousins system (magnitudes B, V, R, T): 

B-b = 26.14-b 0.16(6-u) 

V-v = 26.30- 0.11(6-u) (1) 

R—r = 26.62 — 0.02(u — r) 

I-i = 25.77-b0.07(r-f) . 

The instrumental magnitudes were calculated as 

6 = -2.51og(F{,/texp) - AaTOfc - kb/ cosZ , (2) 

etc, where F is the DN value in a given aperture, texp 
is the exposure time, Aam is the correction for the finite 
aperture (obtained from the Point Spread Function [PSF] 
of bright stars), and k/ cosZ is the correction for the at¬ 
mospheric extinction. Errors of the numerical coefficients 
in eq. (1) do not exceed ±0.01. 

The errors of the stellar magnitudes can be determined 
by measuring the signal-to-noise (S/N) ratios: 

Am = -2.5 log(l ± (3) 


1 through 11 in Fig. 1), together with the effective wave¬ 
lengths and widths of the Cousins filters, are presented in 
Table 3. Figure 2 shows the comparison of the color in¬ 
dices for 39 objects of the pulsar field, from V = 17.5 to 
V = 22.5, with the color indices for stars from the Lan¬ 
dolt catalog. The color-color dependences are close to each 
other, which attests both low absorption in the Calaxy for 
the direction towards PSR B0656±14 and the absence of 
systematic shifts in our photometric system. 

2.5. Optical counterpart of PSR B0656+14 

Contour plots of the pulsar vicinity in all the 4 filters ob¬ 
tained from smoothed (3x3 pixels) images are shown in 
Figs. 3 and 4 The first five isophotes (for lower bright¬ 
nesses) correspond to the following values above the back¬ 
ground: Sky + n\JSky/Npix, with n = 1, 2, 3, 4, 5, where 
Npix = 28 corresponds to the aperture diameter of 6 pix¬ 
els (1'.'6), and the value of Sky is that for the background 
near the pulsar. The isophotes for higher brightnesses are 
scaled arbitrarily to demonstrate brighter details. It is seen 
from Figs. 3 and 4 that 5 objects are detected within < 5" 
around the expected pulsar position. The position of one 
of them (“oljPSR” in Table 2 and Fig. 6 ) coincides, within 
the budget of astrometry errors, with the radio position, 
so that this object is the most plausible candidate for the 
optical counterpart of PSR B0656±14. The nearest object 
“o2” (detected in the I filter only) lies to the North at a 
distance of « 2". Since the objects in the near vicinity 
are not seen simultaneously in all 4 filters, we have shown 
them, for the sake of convenience, in the separate scheme 
of Fig. 6 . 

The photometry of the pulsar candidate (Table 4) was 
carried out for several aperture radii Ra, to optimize the 
S/N ratio. The object “o2”, whose brightness in the I fil¬ 
ter is comparable to that of the pulsar candidate, starts to 
influence the photometry of the candidate at Ra > 4 pix¬ 
els ( > 1"). Most notations of Table 4 have been described 
above. The columns mins and mcous are for the magnitudes 
in the instrumental and Cousins systems, respectively. The 
last column shows the source flux (in 10 “^® erg cm“^ s“^ 
Hz“^) obtained from the Cousins magnitudes with the aid 
of the data on aLyr by Fukigita et al. (1995) 


Fv Gain 
^F + Sky X Np 


( 4 ) 2.6. Objects in near vicinity of the PSR B0656+14 coun¬ 

terpart 


where Npi^ is the number of pixels in the aperture chosen, 
Sky is the value of the local background in DN per pixel, 
and Gain = 2.3e“/DN is the gain. The CCD readout 
noise is ~ 8 e“ or 4 — 5 DN and one was negligible in our 
observations. 

Making use of eqs. (l)-(4), we measured the magni¬ 
tudes of relatively bright stars of the pulsar field (indi¬ 
vidual images were summed for each filter). The Cousins 
magnitudes of 11 of these stars (denoted by the numbers 


The coordinates of the 5 objects within < 5" of the ex¬ 
pected pulsar position (see Figs. 5 and 6 ) are presented in 
Table 2. 

The object “ol,PSR”, whose brightness distribution is 
consistent with that of a point source, is the only probable 
candidate for the pulsar counterpart, based upon both its 
position and photometry. 

The object “o3” is seen in all our B, V,R,I images, as 
well as in the contour plots of Fig. 2 of Caraveo et al. 
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Table 3. Stellar magnitudes of 11 stars around the pulsar. 


BVRI 
(Aefi = 4448 A, (Aeff = 5505 A, (Aeft = 6588 A, (Aes = 8060 A, 
FWHM=1008 A) FWHM=827 A) FWHM=1568 A) FWHM=1542 A) 


1 

19.35 

± 

0.01 

18.31 

± 

0.01 

17.72 

± 

0.01 

17.17 

± 

0.01 

2 

18.71 

± 

0.01 

18.02 

± 

0.01 

17.66 

± 

0.01 

17.32 

± 

0.01 

3 

19.90 

± 

0.01 

18.69 

± 

0.01 

17.96 

± 

0.01 

17.36 

± 

0.01 

4 

20.94 

± 

0.02 

19.84 

± 

0.01 

19.19 

± 

0.01 

18.64 

± 

0.01 

5 

22.61 

± 

0.02 

21.31 

± 

0.01 

20.54 

± 

0.01 

19.93 

± 

0.01 

6 

22.70 

± 

0.03 

21.96 

± 

0.02 

21.52 

± 

0.01 

21.12 

± 

0.02 

7 

21.99 

± 

0.03 

20.31 

± 

0.02 

19.13 

± 

0.01 

17.74 

± 

0.02 

8 

23.08 

± 

0.03 

21.52 

± 

0.02 

20.50 

± 

0.01 

19.53 

± 

0.02 

9 

23.58 

± 

0.03 

22.06 

± 

0.02 

21.24 

± 

0.01 

20.44 

± 

0.02 

10 

24.52 

± 

0.03 

22.88 

± 

0.02 

21.92 

± 

0.01 

21.08 

± 

0.02 

11 

21.38 

± 

0.03 

20.40 

± 

0.02 

19.89 

± 

0.01 

19.39 

± 

0.02 


<— 39 Objects From P5R0656+14 Field 
<— stars From Landolt's catologue 


. + 4 ^ 


< - 39 objects From PSR0656+14 Field 

< - stars From Landolt's catologue 


+ 4 +? 

■ ■ ,fA+- 

0.5 


Fig. 2. a) V — R vs. B — V diagram for stars of the PSR B0656+14 field. Crosses denote the color indices for the field objects. 
Points show the color indices for stars from the Landolt catalog, b) R — I vs. V — R diagram for the field stars. 


(1994a). It has been also detected as a faint extended 
object in the HST F130LP image (Pavlov et al. 1996a). 
Moreover, it can be seen even in our first B images (Kurt 
et al. 1997a) as a background enhancement exactly at the 
same position. Thus, the presence of this faint (extended?) 
object at ~ 4" to the North of the pulsar is firmly estab¬ 
lished. The object “o4” has also been detected with the 
ESO-3.6m and NTT telescopes in the V filter (Caraveo et 
al. 1994a), although it is difficult to distinguish between 
the objects “o4” and “o5” because of their obviously ex¬ 
tended nature. These objects are outside of the 7?4 x 7'.'4 
frame of the JIST FOG image presented by Pavlov et 
al. (1996a). 


The object “o2”, closest to the pulsar, is reliably de¬ 
tected only in the / filter. In other filters only a non¬ 
significant enhancement of the background is seen at this 
position. We suppose that the object “o2” is not seen in 


the HST FOG image because it is too red to be detected 
in the F130LP band. 


Table 5. Stellar magnitudes of objects around 
PSR B0656-bl4. 


object 

B 

V 

R 

I 

o2 

- 

- 

- 

23.83 

(±0.2 1;-0.19) 

o3 

25.0 

±0.2 

25.4 

±0.4 

24.63 

(±0.15;-0.13) 

23.93 

(±0.2 4;-0.20) 

o4 

25.5 

±0.35 

25.5 

±0.5 

24.59 

(±0.21;-0.17) 

23.50 

(±0.12 ;-0.10) 

o5 

- 

- 

24.35 

(±0.14;-0.12) 

23.54 

(±0.20;- 0.17) 




















Fig. 5. Vicinity of the pulsar PSR B0656+14 in the B, V, R, I filters. The size of each fragment is 54" x 54". The optical 
counterpart of the pulsar (“of,PSR”, cf. Fig.6) is indicated by arrows in the B and / images. The nearest object (“o2” in Fig. 6) 
lies to the North at a distance of ss 2". 


Table 5 gives estimates for magnitudes of these objects. 
These estimates and the distances of the objects relative 
to “of,PSR” (Table 2) allow us to conclude that these 
objects do not affect the measured values of the B, V, 
R magnitudes of the pulsar counterpart. A minor effect 
about 3% from whole flux of the nearest object “o2” on 
the I magnitude determined for the pulsar counterpart 


with the aperture radius of 3 pixels (diameter of 1V6) can 
be expected, but it should not exceed the errors of our 
measurement. 
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Table 4. Photometry of the PSR B0656+14 optical counterpart. 


filter 

Sky /texp 

(DN/pix/s) 

Ra 

(pix) 

Npix 

S/N 

-2.5l0g(F/texp) 

Aain 

k/ cos Z 

^9 

^^COUS 

Flux 

(fiJy) 



3 

28 

5.6 

-t0.13 

1.01 





B 

2.465 

4 

50 

6.1 

-0.24 

0.63 

0.43 

-1.30 

24.85(+0.19,-0.16) 

0.46 ± 0.08 



5 

79 

6.1 

-0.42 

0.41 







3 

28 

6.6 

-0.29 






V 

6.429 

4 

50 

6.7 

-0.58 

19 

0.26 

-1.35 

24.90(+0.16,-0.14) 

0.40 ± 0.06 



5 

79 

7.1 

-0.73 








3 

28 

7.9 

-1.11 

ESI 





R 

18.207 

4 

50 

9.8 

-1.34 

19 

0.19 

-2.09 

24.52(+0.12,-0.11) 

0.47 ±0.06 



5 

79 

9.3 

-1.60 

IBI 







3 

28 

4.6 

-1.02 

0.81 





I 

29.922 

4 

50 

4.2 

-1.42 

0.49 

0.12 

-1.95 

23.81(+0.27,-0.21) 

0.71 ±0.16 



5 

79 

4.1 

-1.67 

0.31 







Fig. 6. Objects from nearest vicinity of the pulsar are indicated 
whose coordinates are given in Table 2, and estimates of stellar 
magnitudes are given in Table 5. Coordinates (in arcseconds) 
correspond to those in Figs. 3 and 4. 


3. Discussion. 

Figure 7 shows the UV-optical energy fluxes in different 
bands for the PSR B0656+14 counterpart obtained with 
the HST and 6 -m telescope. It demonstrates that our re¬ 
sults (crosses B, V, R and I), being compatible with pre¬ 
viously obtained data, considerably extend redward the 
observed frequency domain, fill in the gap between the 
visual and UV parts of the spectrum, and provide more 
accurate fluxes in the V and B bands. The overall spec¬ 
trum firmly confirms a nonthermal origin of the optical 


radiation suggested by Pavlov et al. (1996a). On the other 
hand, the relatively high value of the F130LP flux makes 
very plausible that the thermal radiation from the NS sur¬ 
face contributes noticeably into the UV part of the spec¬ 
trum. Parameters of the thermal and nonthermal fluxes 
can be constrained by fitting the observed spectrum with 
a two-component model (cf. Pavlov et al. 1996a), a sum 
of the power-law and thermal spectra (the latter has the 
form of a Rayleigh-Jeans spectrum at these frequencies for 
plausible NS temperatures): 


fH 





(5) 


where t'o is an arbitrary reference frequency (we chose 
t'o = 8.766 X lO^"* Hz, which corresponds to Aq = 3420 A), 
A{i/) is the interstellar extinction (see, e. g., Savage & 
Mathis 1979), /o and go are the values of the nonthermal 
and thermal fluxes at v = i>o- The latter can be expressed, 
for the chosen value of i^o, as 


go = 3.116 X 10"^^G 


erg 

cm^ s Hz ’ 



where T = lO^Te K is the NS brightness temperature, 
i?oo = 10i?io km is the NS radius (as seen by a dis¬ 
tant observer), and d = SOOdsoo pc is the distance to 
the pulsar. According to Pavlov et al. (1996a), a plausi¬ 
ble value for the color excess, which determines the ex¬ 
tinction curve A{v), is E{B — V) = 0.03 {Ay = 0.09). 
We fixed E{B — V) at this level and fitted the fluxes 
in the 6 bands {B, V, R, I, F555W and F130LP) with 
eq. (5), varying /q, a and G. The best-fit model (min¬ 
imum = 4.69 for 3 degrees of freedom) corresponds 
to fo = 1-75 X 10“^° erg cm“^ s“^ Hz“^, a = 1.53 and 
G = 4.1. The best-fit spectrum is shown in Fig. 7 with the 
solid line, and the best-fit thermal and power-law compo¬ 
nents with dashed lines. 



























































10 


V. G. Kurt et al.: BVRI observations of PSR B0656+14 



Log Frequency (Hz) 

Fig. 7. Fluxes observed from the PSR B0656+14 optical coun¬ 
terpart in different spectral bands. The crosses B, V, R and / 
show our results; the crosses 130LP and 555W show the fluxes 
detected with the HST FOC/F130LP by Pavlov et al. (1996a) 
and WFPC2/F555W by Mignani et al. (1997). Dashed crosses 
show the fluxes obtained by Caraveo et al. (1994; V band) and 
Kurt et al. (1997a; B band). The 3a upper limits on the flux 
from undetected sources in the B, V,R,I filters are shown for 
the aperture diameter of 6 pixels (« l”6). The solid line shows 
the best fit of the fluxes in the 6 bands with the two-component 
model (power law -|- Rayleigh-Jeans spectrum; eq. [5]) for the 
color excess E{B — V) = 0.03. The dashed lines are the com¬ 
ponents of this model. 

The value of G can be compared with what is ex¬ 
pected from extrapolations of the soft X-ray spectral flux 
observed with the ROSAT and ASCA (Finley, Ogelman 
& Kiziloglu 1992; Anderson et al. 1993; Greiveldinger et 
al. 1997) to the optical range. The result of this extrap¬ 
olation depends on the model assumed to fit the X-ray 
spectrum. The blackbody fit of the ROSAT spectrum 
yields Tq = 0.9 ± 0.04 for i?io/d 5 oo = 1-0 ± 0.2 (Finley 
et al. 1992), which corresponds to the Rayleigh-Jeans pa¬ 
rameter G = 0.9 ± 0.3. A fit of the same data with a set 
of hydrogen NS atmosphere models (Pavlov et al. 1995) 
gives, for B = 4.7 x 10^^ G, the radius-to-distance ra¬ 
tio Rio/dsoo = 2.3 ± 0.4 and the effective temperature 
= (5.3 ± 0.5) X lO^K (Anderson et al. 1993). Since 
the model spectrum differs from the blackbody spectrum, 
there is no a universal relation between the effective and 
brightness temperatures; however, the model spectrum in 



Fig. 8. Confidence levels (20, 68, 90 and 99%) for fitting of the 
observed spectrum of PSR B0656-I-14 with the two-component 
model (eq. [5]) 

the optical range closely reminds, and can be fitted with, 
the Rayleigh-Jeans spectrum with Tq ~ 0.9T^, thus giv¬ 
ing G ~ 2.6 ± 0.7, almost thrice greater than G inferred 
from the blackbody model. 

We see that the best-fit value of G obtained from the 
available optical-UV observations exceeds the predictions 
of both the blackbody and hydrogen atmosphere interpre¬ 
tations of the soft X-ray flux. However, as Fig. 8 shows, 
the fitting parameter G is constrained poorly by our data: 
0<G<6.2atal(T level. The main reason of that is the 
lack of data in the far-UV range where the thermal com¬ 
ponent gives the main contribution. Thus, we can only 
conclude that the results are consistent with the presence 
of a thermal component, but we cannot distinguish be¬ 
tween different NS surface models and can only estimate 
upper limits for the Rayleigh-Jeans parameter, G < 6.2, 
6.8 and 7.3 at the confidence levels of 68, 90 and 99%. 
We expect that better constraints will be obtained from 
recent HST FOG observations of this pulsar in the UV 
range (Pavlov, Welty & Cordova 1997). 

Figure 7 clearly demonstrates that the pulsar radia¬ 
tion in the IR-optical range is predominantly of a nonther- 
mal origin, in agreement with the conclusion of Pavlov et 
al. (1996a). An independent confirmation of the nonther- 
mal nature comes from the high pulsed fraction observed 
by Shearer et al. (1997) in the B band. Although the slope 
of the nonthermal spectrum is not well constrained be¬ 
cause of relatively large errors (especially in the I band): 
0.3 < a < 2.6 at a 1 (J level, it looks much steeper than 
those of the (much younger) Crab pulsar, a = —0.11±0.13 
(Percival et al. 1993; see also Ransom et al. 1994). The 
nonthermal component of the younger Vela pulsar has the 
slope a ~ 0.3 in the VBU range, but its R flux is about 
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twice lower than the extrapolation of the power law (Na- 
suti et al. 1997). Even more drastic difference is seen be¬ 
tween the shapes of the optical spectra of PSR B0656-I-14 
and the older pulsar Geminga (r = 3 x 10® yr) which 
looks very much alike B0656-I-14 in the soft X-ray range. 
Geminga’s flux has a minimum at the B band, is peaked 
at the V band, and falls off sharply redward, which al¬ 
lowed Bignami et al. (1996) to suggest the presence of 
a proton cyclotron line centered at A ~ 6000 A. On the 
contrary, the flux of PSR B0656-I-14 grows redward of V, 
thus demonstrating that the properties of the nonthermal 
IR-optical-UV radiation are neither the same in different 
pulsars nor they have a clear correlation with the pulsar 
age. 

The apparent excess of the B flux with respect to ac¬ 
ceptable continuum fits (Fig. 7), responsible for the rela¬ 
tively large value of the minimum may indicate the 
presense of a spectral feature at these wavelengths, per¬ 
haps similar to that reported for Geminga. More observa¬ 
tions are needed to verify whether the feature is real. But 
at least three broad emission lines above 3600 A were re¬ 
cently detected indeed in a prism spectrum of PSR B0540- 
69 with the FOS/HST by Hill et al. (1997). The authors 
connect directly the presence of these strong emissions 
with a nebula surrounding B0540-69. 

The nonthermal component has also been observed in 
high-energy radiation from PSR B0656-I-14. Greiveldinger 
et al. (1997) reported a power-law component in the X-ray 
spectrum of PSR B0656-I-14, with a = 0.5 ± 1.1, which 
dominates at if > 2 keV. Ramanamurthy et al. (1996) pre¬ 
sented an evidence for y-ray emission (E > 100 MeV), 
with a = 1.8±0.3. Being extrapolated to the X-ray range, 
our best-fit power-law component {a = 1.53) gives the X- 
ray flux about 3 orders of magnitude lower than observed 
by Greiveldinger et ah, i. e., the nonthermal flux should 
grow somewhere between the UV and X-ray ranges. To 
merge with the X-ray power-law component at if ~ 2 keV, 
our optical flux should have a slope a ~ 0.3 — 0.4, hardly 
compatible with our fit (see Fig. 8), and in this case the 
spectrum should steepen at higher energies in order to be 
compatible with the y-ray data. It is, however, hard to 
believe that the spectrum would have the same slope in 
the optical and high-energy ranges — in fact, none of the 
four other pulsars with investigated nonthermal optical 
components shows such behavior. 

To put more stringent constraints on the parameters 
of the thermal and nonthermal components, more obser¬ 
vations are needed. Observations in the IR range would be 
most useful for elucidating the slope a of the nonthermal 
component, whereas far-UV observations would be crucial 
to measure the Rayleigh-Jeans parameter G proportional 
to the temperature of the NS surface. Since G and a are 
correlated with each other (see Fig. 8), measuring one of 
these parameters with a higher accuracy would immedi¬ 
ately constrain the other parameter. It would be also use¬ 
ful to obtain more detailed spectrum around the B band 


in order to understand whether the excess we observed 
is associated with a spectral line. Such observations are 
quite feasible with both the HST and large ground-based 
telescopes. 
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